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^ ^ I Abstract. Corrections to solar system gravity are derived for f{Q) gravity theories, 

in which a function of the Gauss-Bonnet curvature term is added to the gravitational 
action. Their effects on Newton's law, as felt by the planets, and on the frequency shift 
of signals from the Cassini spacecraft, are both determined. Despite the fact that the 
Gauss-Bonnet term is quadratic in curvature, the resulting constraints are substantial. 
It is shown that they practically rule out f{G) as a natural explanation for the late-time 
acceleration of the universe. Possible exceptions are when f{G) reduces to something 
1^ . very close to a cosmological constant, or if the form of the function / is exceptionally 

fine-tuned. 
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^ ' 1. Introduction 

^ ■ The current accelerated expansion of our universe cannot be explained by conventional 
^ ■ general relativity if our universe contains only standard matter and radiation. Some 
form of additional dark energy, such as a cosmological constant, may be the source of 



the acceleration, although such models suffer from serious fine-tuning problems. It may 
be that the acceleration instead comes from corrections to Einstein gravity. Such an 
approach has the potential to avoid the fine-tuning problems, resulting in a far more 
credible theory. One candidate for effective dark energy is the quadratic curvature 
Gauss-Bonnet term 

g = R^- AR^.R^'' + R^^p^R^""'' , (1) 

which is a natural addition to the Einstein- Hilbert action [1]. 

In fact, on its own in four dimensions, the contribution of the Gauss-Bonnet term to 
the gravitational field equation is trivial. For it to have an effect the theory needs to have 
extra dimensions, such as in the brane world scenario with one [2\ of more [3] additional 
dimensions, or alternatively the Gauss-Bonnet term can be coupled to a scalar field. 
Another related possibility is to add a function /(^) of it to the gravitational action [1]. 
Such a theory will be the subject of this article. 

Its potential to give a more elegant explanation for the universe's accelerated 
expansion makes gravity modification very appealing. However it does have one major 
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drawback, namely that deviations from general relativity will be felt on all scales, not 
just cosmological ones. In particular gravity will be altered within our solar system, 
where high precision tests of general relativity have been performed. As we will show 
for f{Q) gravity, the corrections are generically too large, and allow it to be ruled out 
as a solution to the dark energy problem (with some severely fine-tuned exceptions). 

Within our solar system, gravitational fields are weak. Since the Gauss-Bonnet term 
is quadratic in curvature, it might be expected that its effects will be sub-dominant, 
and not significantly constrained. However as we will demonstrate, this reasoning is 
flawed. In section [2] we will show, after introducing the theory, that if the cosmological 
contribution of f{Q) is to be large enough to act as dark energy, the couplings in / 
must be extremely large. This will greatly magnify the effects of Q in the solar system, 
and produce corrections to the Newtonian and post-Newtonian potentials, which we 
will derive in section [3l Strong constraints on these corrections arise from planetary 
motion and light bending measurements, which are analysed for a more general theory 
in section m In section [5] these are applied to f{Q) gravity, and it is shown that the 
gravitational effects within the solar system are indeed large enough to conflict with the 
observational data that is available, and thus allow this dark energy candidate to be 
practically ruled out. 



2. Dark energy from f{Q) gravity 

Working in units with c = 1, we will be studying the theory 



(2) 



It is equivalent to coupling a Gauss-Bonnet term to a scalar field with a potential, 
but without an explicit kinetic term (in contrast to the Gauss-Bonnet theories studied 
in [5], [6]). In particular, if we take the action 

C = ^^[R + a<P)G-2V{<P)] (3) 

with 

y(0) = -/(0) + 0/(0) , = 2/'(0) , (4) 

and then vary it with respect to the scalar field, we obtain (p = Q. Substituting this 
back into the action ([3]), it reduces to the theory ([2]). 

Varying the action ([2]) with respect to g^^, we obtain the gravitational field 
equations 

+ 9,{Rppya + 2Rp[i,g^]n — 2Rp[^g„]p + Rgp\^,g„^p)V^V f {Q) 

These can also be obtained from ([2]). The constant Gq in the above equation is the 
gravitational coupling of matter, which may or may not be the same as the gravitational 
coupling G that we perceive on Earth. 
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In a cosmological background (rfs^ = —dt'^ + a{tydxl), we find Q = 2AH'^ {H + H'^) , 
with H = a/a. The Friedmann equation can be written as 

= ^^Pmat + i^g , (6) 
where the Gauss-Bonnet density fraction is 

% = -8HdJ'ig) + ^I^^^. (7) 
We will focus on the case f{Q) = C^". The expression then reduces to 

1 + 3w 

where w = —1 — 2H/{3H'^) is the effective equation of state for the universe. An 
important point to note is that if fig ~ 0.7, as is required if f{Q) is to give a sufficient 
contribution to the dark energy density, then C ~ Hq~^^. This is extremely large (for 
n > 1/2), and so we see that the Gauss-Bonnet term must be very strongly coupled if 
it is to have any chance of explaining the accelerated expansion of our universe. 

Note that if a solution of the above theory is solve the dark energy problem, 
Qg ~ 0.7 is necessary, but not sufficient. We also need the solution to actually produce 
enough acceleration (so w ~ —1), and for the cosmological evolution of our universe to 
reach it (after passing though a period of matter domination, exactly like the one that 
occurred in our universe). In this work we will mainly be concerned with the magnitude 
of fig, which will be enough to rule out most possible f{Q) dark energy models. 



3. Solar system gravity 

Within the solar system we can use the approximate Post-Newtonian metric [7] 

ds^ = -(1 + 2$) dt^ + (1 - 2^) Sijdx'dx^ + 0(6^/2) _ (g) 

with $, ^ ~ e and $, ^ ~ e^/^. The gravitational fields of the sun and planets are 
relatively weak and slowly varying, and e < 10~^. This allows us to study gravity with 
a perturbative expansion in e. Note that the metric is motivated by the properties 
of our solar system, and not by the choice of gravitational theory. 

Making no assumptions about the size of / or its derivatives, we find that to leading 
order in e, the gravitational field equations are 

A^ = AnGop^.t+fiG)-QnQ)-mfiG),^+'^)+^'oii,gf,g'n+eOif,gniio) 

AM/ = 4nGoPm.t + ^^^^ - AV{f\g), + 6^0(1, gf\ g'f") + e 0(/, ^/')(11) 

and g = 8X'($, \I') + O(e^). We have introduced the operators 

AF = J2 F,^^ , '^(X, Y) = Y, X,^,Y,, - AX AY , (12) 

i id 
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which for functions with only r dependence reduce to 

AF = \dr{r'^drF) , Y) = -^dr{rdrXdrY) . (13) 

In principle, the higher order (in e) terms of the above expansion could give significant 
contributions if their coefficients are very large. However for / = CQ", the corresponding 
coefficients are all of comparable size, and the higher terms can be dropped. Furthermore 
when / = CQ"', we see that the f{G) and Qf'{Q) terms are higher order in e than the 
T>{- ■ ■) terms. For the rest of this section we will assume this form for the theory, and 
so can drop those two terms from the field equations. In this work we will be mainly 
interested in > 0, since theories with negative n have been studied elsewhere in the 
literature. In particular, it has been shown while they satisfy solar system constraints [8], 
they are unable to solve the dark energy problem [9]. We will discuss these points in 
more detail in section [5l 

The usual, Einstein gravity (/ = 0) solution is $ = \1' = — t/, with 



U = 47rGo / rf^x^ ™^ J/ . (14) 







\x - 


- x'\ 



Treating the sun as a uniform sphere we find that t/, and the Gauss-Bonnet term, reduce 
to 

t4xt = — - — , ^cxt ~ 48 — (15) 

outside the sun (r > -R©), and 

Gome/ r2 \ (GprnQ)^ 
= ^ I' - ^ J ' " "''^^ ^''^ 

inside it (r < -R0). 

A general analysis for non-trivial / would be difficult. However we know that 
for any viable theory, the resulting gravitational potentials must be very close to the 
standard form. Following the approach of {5], we start by assuming that the potentials 
are very close to the usual 1/r form. Deviations from this can be treated perturbatively. 
By bounding their size, we can derive constraints on /. Of course for a wide range of 
f{Q) theories, $ and \l/ will be radically different from 1/r, and our approach will not 
give their approximate form correctly. However, since we already know (by definition) 
that such theories fail to give the observed gravitational potentials, and there is no need 
to study them in the first place. Hence our approach will cover all potentially relevant 
cases. 

For the interior of the sun, we see that I^(/'(^int), ^int) = 0. The interior 
solution ([ISD, up to the addition of a constant, is therefore still valid (recall that we are 
ignoring the /, Qf terms as they are sub-dominant for / = CQ"'). 

For the exterior solution we take $ = —CmQ/r + 6^ and \E' = —GniQ/r + S'^. Here 
G is the approximate gravitational coupling that we perceive, and need not be equal to 
the bare coupling Gq. The perturbations to the usual potentials resulting from the / 
dependent T>{- ■ ■) operator are given by 5$ and 6"$. Such an analysis will only be valid 
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if 5$, 5^ <^ U, which will indeed be the case for us. The field equations (fTU]) and (fTTi) 
reduce to 



A5<^^8V{ CuQ 



->n-l 



Gnirr 



Q 



A 5^ ^ 41) Cng 



->n-l 



© 



■ ^^^^ 

with Q ^ 48(Gm0)^/r^. Solving these, and requiring continuity of $, and their first 
derivatives, gives the perturbed exterior potentials 



r 
r 



Gi 



GA 



1 + s) 
[Gi-GAr\nr] 



_m0 
r 

r 



GB 
— ^ — 

(l + s)r« 

[Gi - GBr\nr] 



where we have introduced 

A = 2B = 2n{n - 1) GAS^'iGmQ^^''-^^ 



G\ — Gq + ARr: 



' 



(for s ^ -1) 

(for s = -1) (18) 

s = 2(3n-2).(19) 



We see that there are finite width effects appearing in the gravitational potential. This 
is due to the presence of higher than second order derivative operators appearing in the 
gravitational field equations. This also implies that the usual treatment of the sun, and 
other objects, as point-like masses is no longer valid. 

Clearly the above expressions f|T8|) are different from the usual Einstein gravity 
result, and are likely to come into conflict with the high precision measurements of 
gravity in the solar system. By bounding the corrections, we will be able to constrain 
the parameters in the function f{Q), and hence its dark energy contribution. 



4. Solar system gravity tests 

We will now review two sources of gravitational constraints coming from the solar 
system. First a test of Newton's law from planetary motion, and then the frequency 
shift of light rays, which tests a relativistic effect. The results will apply to any theory 
giving potentials of the form (ITSil . and not just f{Q) = GQ^ gravity. We will apply 
them to the f{Q) dark energy models in section [5l 

Corrections to the Newtonian potential can be bounded by considering their effect 
on bodies orbiting the sun. For an expression of the form (fTSjl . the gravitational 
acceleration experienced by a body at distance r from the sun is 



dr 



Gi-A 



G 



Gmcsir) 



(20) 



For a body following an elliptical orbit with semi-major axis a, Kepler's third law gives 
the period of the orbit as 2t:^ j (GniQ). Bounds on corrections to the effective mass 
are then related to the uncertainties in the measurement of a by 



(21) 







where the index a runs over all bodies which are orbiting the sun. Values of Sa for the 
planets can be found in [10]. The above relation has previously been used to bound 
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dark matter in the solar system [TT], the cosmological constant [12], and another class 
of Einstein-Gauss-Bonnet gravity models [5]. 

As it stands, the above constraint ( l2Ti) depends on Gi, which is related to the 
unmeasured constant Gq. By combining constraints from two bodies, we can eliminate 
it to obtain 

|^(a„--a-)|<3f^ + ^V (22) 



da dp 

which needs to be satisfied for all choices of a, (3. For s > — 1 the strongest bounds 
come from the inner planets, firstly because these are the bodies for which we have the 
best data, and secondly since they are closest to the sun, where the gravity corrections 
are strongest. 

To satisfy the bound fl22|) . either A will need to be very small, s will need to be close 
to zero (in which case the corrections to the potential will mimic Newtonian gravity), 
or s will need to be less that —1 (in which case the corrections are suppressed at the 
smaller distances where better data is available). 

Further constraints come from signals between man-made spacecraft and the Earth. 
The sun's gravitational field produces a time delay in the signals, measurement of which 
provides an additional test of gravity in the solar system. Furthermore, unlike the 
planetary constraint (122|) . this is sensitive to relativistic effects, and so can be used to 
rule out models which mimic Newtonian gravity. 

For a light-ray starting at the Earth, passing close to the sun's surface, continuing 
to the spacecraft {r^ from the sun), and then returning by the same route, the time 
delay in the signal is 

At = -2 [^{r) + ^{r)]dx. (23) 



The signal's path is approximated by r = \^aP+lP, where the impact parameter b 
is defined as the smallest value of r on the light ray's path. A small value of b will 
maximise the above time delay. Particularly good data was obtained for the Cassini 
spacecraft while making its journey to Saturn. During 2002, the impact parameter fell 
as low as 6 = I.GRq ~ 0.0074 AU (the value of re at this point was 8.43 AU). The actual 
measurements obtained were not of the time delay, but its frequency shift given by 



dAt dAt db AGrup-, db r , , , ^ 

y = ^- — -— = 2 + 2.1 ± 2.3 X 10"^ . 24 

^ dt db dt b dt^ ^ ' ^ ^ ' 

For potentials of the form ( ITSl) . the frequency shift evaluates to 

y = (2^1X0 -[A + B]GI,) , (25) 
where we have defined 



2 J^^ (62 + a;2){3+^)/2 

which evaluates to a hypergeometric function. A similar analysis was applied to 
a different Einstein-Gauss-Bonnet theory in [5], with s = 6 due to the simplifying 
assumptions imposed on the model. 
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For small impact parameter, h <C Oe^^e, which is the case for the above Cassini 
bound fl2^ . the leading order behaviour of the above integral is 

/^r([2 + s]/2) 1 



I. 



2r([3 + s]/2) 



s > 



(27) 



s <-2 



2(2 + s) 

In particular Xq ~ 1. 

Comparing (!2ll) and (125!) gives a bound on A 
since like fl2Up it depends on the undetermined quantity Gi. To eliminate that, we 
combine the constraint with fl2Tl) . giving 



B, although this is of limited use. 



(2 X 10" 



6- 



6an 



< 2Aa-' - (A + B)Is < (4.4 x 10~^) + 6 



5ar 



(28) 



Note that the above constraint will need to be satisfied for all choices of a. 

Another well-known solar system gravity test comes from the perihelion precession 
of Mercury. However we will not consider it here as the linearised analysis of section [3] 
is insufficient to determine it, and a higher order (in e) expansion is needed. As it 
turns out, the above results will be sufficient to eliminate f{Q) as a viable dark energy 
component. 



5. Constraining f{G) gravity 

Combining the expression for the effective dark energy fraction ([8]) with the planetary 
motion constraint (1221) . and taking the current effective equation of state for our universe 
to be w = —1, we obtain the bound 



\nc\ < 



5a, 



a _^ Sa^ 



ap 



where Hn 



8.8 X 10 AU is the Hubble distance, and r„ 



Gmr. 



(29) 



9.7 X 10-9 AU 



is the gravitation radius of the sun. For n > 1 we find the strongest constraint comes 
from taking (a,/5) to be Mercury (a ^ 0.39 AU, 5a/ a 1.8 x 10"^^) and the Earth 
(a ^ 1 AU, 5a/ a ^ 0.98 x lO'^^)^ jhis implies 

< 1.6 X 10-^^ (30) 

which is clearly far too small to account for the dark energy of our universe. For smaller 
n, a significant dark energy fraction may be possible. Using all the planets in our solar 
system, we find fig < 0.1 unless 



2 

n 

3 



< 1.1 X 10 



'27 



or 



n < 0.074 . 



(31) 



The two numerical values come respectively from the Earth- Mercury and Earth- Uranus 
combinations (Uranus has a ~ 19 AU, 5a/ a ~ 1.3 x 10"^). Note that for the above 
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bounds on Vtg we assumed the equation of state w = — 1 for our universe. Other values 
will give slightly weaker or stronger constraints, but they will be of the same order of 
magnitude. 

Although the first of the exceptions (!3T|) mimics the correct Newtonian limit 
of gravity, it fails to produce the correct relativistic effects. Using the Cassini 
constraint fl28l) . we find that for the parameters f|T9|) . 

2 X 10-^ 5a Thuf^ 4.4 x 10'^ 5a^ 
6 ^ < "^f^ " ^^^"-^^ < 6 + ^ ' ^^^^ 

where we have again used ([8]) with w = —1. Taking a to be Mercury, we find that 
< 1.2 X 10~^° for n > 0.66, which includes the first range in (l3Til . thus ruling it out. 
It should be pointed out that the errors 6aa appearing in [10] are in fact the 

statistical errors in the values of coming from a least-squares fit of observations. 

The real errors may well be an order or two of magnitude higher, and would give 

correspondingly weaker bounds. However the above bounds are so strong that even this 

would not be enough to save the theory. 

The second exception (!3T!) to the bound (129|) includes the parameter range n < 0. 

Such models have been studied elsewhere, for example [8], [9l [H]. We see that (for 

n < 1/2) the correction to the Newtonian potential for a general mass M is 

(3 \ 1— 2?i 

which actually decreases for larger masses. Obviously, the above expression is not valid 
at very large distances, where 5$ will cease to be a small perturbation. However for 
couplings of the size needed to solve the dark energy problem, C ~ ifg""^", the above 
5$ will remain small within the solar system. 

We see that for a heavy object such as the sun the corrections to $ and \1/ are 
suppressed, which is why (for low enough n) they do not conflict with our solar system 
constraints. On the other hand, we see that the above expression becomes large for 
small masses, implying that the weak field approximation we have been using breaks 
down. This suggests that the theory will predict non-Newtonian gravity for table-top 
gravity experiments, and so will disagree with observation. However it must not be 
forgotten that laboratory experiments are performed on the Earth, which has a large 
gravitational field of its own. The above expression ( l33l) no longer applies since it is 
now the Earth which gives the dominant contribution to the Gauss-Bonnet term Q. The 
resulting bounds on the theory are then far weaker, and do not rule out dark energy 
models with n < P [Hj. 

Although f{Q) gravity with inverse powers of Q can give a large enough Qg, and 
avoids conflict with solar system gravity tests, this is not enough to solve the dark energy 
problem. During its early evolution, our universe passed through a matter dominated, 
decelerating phase (with Q = 2AH^{d/a) < 0), before entering the current accelerated 
phase (with ^ > 0). The Gauss-Bonnet term must therefore have passed though zero 
at some point. We see from the cosmological field equations ([6]), ([7]) that if /(O) is not 
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bounded, ^ = corresponds to curvature singularity, and is unreachable (for finite H). 
Since no transition to accelerated expansion occurs, inverse powers of Q do not give 
viable dark energy models [9]. However, they may still be relevant to the dark matter 
problem |14j . 

So far we have only considered a special subset of the modified gravity theory ([2]) 
where the function f{Q) is some power of the Gauss-Bonnet term. A more general 
theory can be expanded as a power series 



where the n are not necessarily integers. Now each term gives a correction to the 
Newtonian potential $ with a different r dependence, and so barring extreme fine- 
tuning of the Cni we can constrain each of them separately. Similarly for Hence, 
using the above results, we find that for each of their contributions to the dark energy, 
fig"'' <C 1. The only exception is < n < 0.074, which is very similar to a cosmological 
constant. 

Finally, we note that our gravity corrections have turned out to be tiny, at least for 
f{Q) = CQ^ theories with positive n, and which also satisfy the solar system constraints. 
Hence the leading order corrections 5$ and 5'^ used in section [3] are much smaller than 
GniQ/r, and the perturbative analysis used there is indeed self-consistent. In the case of 
= 0, the leading order non-zero corrections to $ and \1/ are not covered by our analysis. 
This case is the same as a cosmological constant, and has been covered elsewhere in the 
literature (e.g. [12])- It is worth pointing out that a very wide range of f{Q) theories 
do not have $ ^ \1/ —GniQ/r as an approximate solution, and so the analysis of this 
paper does not apply to them. However, since the gravitational field of the sun in these 
theories has no resemblance to what we have observed in the solar system, they are 
already ruled out, and there is no need to apply this paper's analysis. 

6. Conclusions 

We have seen that a careful analysis of solar system gravity provides a powerful probe 
of higher curvature gravity theories, and can be used rule out modified gravity models 
which are intended to solve the dark energy problem. Specifically, strong constraints on 
Newton's law are obtained from the motion of the planets. These are sufficient to rule 
out a large fraction of f{Q) dark energy models. One class of model that survives these 
constraints are those with inverse powers of the Gauss-Bonnet term Q, although their 
cosmological evolution is not consistent with our universe. 

The remaining f{Q) models which survive the Newtonian constraints either mimic 
Newtonian gravity, or are very close to a cosmological constant. Further strong 
constraints come from the measurement of the frequency shift of signals form the Cassini 
spacecraft. These rule out the models which mimic Newtonian gravity. Hence f{Q) is 
practically ruled out, as it is only viable when it is virtually identical to a cosmological 
constant. 




(34) 



n 
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In principle, with enough fine-tuning, a model could be constructed where the 
deviations from general relativity cancel for each of the planets (and for the Cassini 
bound). The non-standard gravity effects could then still be large enough to account 
for the dark energy. However for this to work the corrections must cancel for every 
choice of position r, and mass, for which there is a planet, moon, asteroid, laboratory 
experiment, time delay or frequency shift measurement. Even then, there is no guarantee 
that the non-standard cosmological evolution will be accelerating, be free from ghosts [9], 
and satisfy all other cosmological tests, such as producing correct growth of density 
perturbations, which is also a problem for f{Q) gravity [15j. The problem of ghosts 
could even arise within the solar system, giving another way to constrain the theory 
there. The presence of ghosts can be determined by analysing perturbations of the solar 
system metric. However we will not pursue this for our / = CQ^ dark energy (n > 0) 
solutions, since they have already been ruled out. Furthermore, we did not actually 
obtain a solution to the field equations within the solar system (expect when VLg is 
neghgible), so we do not have a background solution whose perturbations are worth 
analysing. 

Following the example of f{R) gravity [16], an f{Q) that reduces to constant for 
large Q and zero at small Q may give a viable dark energy model while satisfying solar 
system constraints. In fact any f{G) whose form gives acceptable corrections for solar 
system curvatures (e.g. a inverse power of ^), but produces acceleration at cosmological 
scales could be viable. However such a theory will require even more fine-tuning than a 
conventional cosmological constant, so it is debatable whether this is an improvement. 

Our strong constraints on f{Q) dark energy rely on the fact that we can directly 
link its local and cosmological properties. This is possible because the behaviour of the 
theory is fully determined by the form of the metric (which is known at both scales) 
and the constant parameters which determine /. As an alternative type of theory we 
could have extended the action ([3]), to include additional kinetic terms for the scalar 
0, as in O E]. These can all have 0-dependent couplings, and can include higher 
order terms as well as the usual (9^0)^. The resulting theory has the same number of 
degrees of freedom as the one studied in this paper, suggesting the same analysis can be 
applied. However, the value of the field is no longer directly determined by the metric: 
(p = Q. Instead it is related to (^^^ by a differential equation, whose solution can include 
integration constants which are determined by cosmological rather than local scales. 
In particular, the value of at the edge of the solar system is (roughly) zero for f{Q) 
gravity, and arbitrary for the more general theories. Hence in the latter case connecting 
the strong local constraints and the cosmological evolution will require further analysis 
of the theory's behaviour on intermediate scales. The theories studied in [5l [6] therefore 
still have the potential to solve the dark energy problem, although of course further 
work is required. Another, more obvious advantage of the above scalar-tensor theories, 
is that they contain more coupling functions than the single one in f{Q). With a greater 
range of theories to choose from, there is more hope of finding one which satisfies both 
local and cosmological constraints. 
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